A Coupled Quantum Otto Cycle 
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We study the 1-d isotropic Heisenberg model of two spin- 1/2 systems as a quantum heat engine. 
The engine undergoes a four-step Otto cycle where the two adiabatic branches involve changing 
the external magnetic field at a fixed value of the coupling constant. We find conditions for the 
engine efficiency to be higher than the uncoupled model; in particular, we find an upper bound 
which is tighter than the Carnot bound. A new domain of parameter values is pointed out which 
was not feasible in the interaction-free model. Locally, each spin seems to effect the flow of heat in 
a direction opposite to the global temperature gradient. This seeming contradiction to the second 
law can be resolved in terms of local effective temperature of the spins. 

PACS numbers: 05.70.Ln, 07.20.Pe 



I. INTRODUCTION 

Quantum generalisations of classical heat cycles have 
now been studied for some years. When the working 
medium is a few-level quantum system, new lines of en- 
quiry open up due to additional features like discreteness 
of stateSjquantum correlations, quantum coherence and 
so on 1-7]. Many models have served to investigate the 
validity of second law of thermodynamics in the quan- 
tum regime [H, The possibility of small scale devices 
and information processing machines [l(| has generated 
further interest into the fundamental limits imposed on 
the heat generation, cooling power and thermal efficien- 
cies achievable with these models [lll - [l3l ] . Quantum ana- 
logues of Carnot cycles, Otto cycles and other brownian 
machines have been analysed fl3. [l5j. Further, both in- 
finite and finite-time (l6l - [23 | thermodynamic cycles 
have attracted attention. 

The quantum Otto cycle which occupies our interest 
here consists of a working substance with hamiltonian 
H and initial density matrix p being manipulated be- 
tween two heat reservoirs (the reservoir temperatures 
satisfy T\ > T 2 ) under two adiabatic and two isochoric 
branches. On the adiabatic branches, the system is as- 
sumed to follow quantum adiabatic theorem and ther- 
modynamic work is defined in terms of the change in 
energy levels at given occupation probabilities. If the 
hamiltonian is changed from H\ to H% by controlling an 
external parameter then the work performed is defined as 
Tr[p(H2 — Hi)]. On the other hand, while traversing the 
isochoric branches, heat is exchanged with the reservoirs. 
Thus if the density matrix of the system changes from pi 
to p2 for a given hamiltonian H, then heat exchanged is 
T^[(p2—pi)H]. As an example, for an effectively two-level 
system whose energy splitting can be varied from Ei to 
Ei , the Otto efficiency has been found to be 1 — -E2/-E1, 
which is bounded from above by Carnot value due to the 
condition E 2 /Ei > T 2 /Ti 0. 

Recently, some authors have studied the role of dif- 
ferent quantum interactions using spin- 1/2 particles in 
a Quantum Otto cycle In particular, the role of 



quantum entanglement has been conjectured using mea- 
sure like concurence and the second law has been shown 
to hold in such models. In this paper, we also investi- 
gate a coupled Otto engine using a 1-d Heisenberg model 
with isotropic exchange interactions between two spin- 
1/2 particles (see Eq. ([1]) below). In 0] the same model 
was analysed, where during the adiabatic steps, the ex- 
change constant J was altered between two chosen val- 
ues (Ji — > J 2 — > Ji), while keeping the external magnetic 
field at a fixed value. From an experimental point of view, 
it is also interesting to investigate a cycle where the ex- 
change constant is fixed and only the magnetic field is 
varied during the adiabatic steps. Further, the uncou- 
pled model cycles considered earlier in literature can be 
taken as a benchmark with which to compare the engine 
performance of the coupled model. 

The paper is organised as follows. In section II, we 
present the quantum model of our working medium, enu- 
merating the energy eigenstates and eigenvalues. In IIA, 
the various stages of the heat cycle are described and 
expressions for heat exchanged with reservoirs and work 
delivered are calculated. It is instructive to develop the 
engine operation based on local description. It is shown 
that all the work is done locally by each spin. In sub- 
sections IIIA and IIIB we develop two cases i) Bi > B2 
and ii) B 2 > B\. The latter case is possible only in 
the presence of interactions. It is observed for this case 
that second law of thermodynamics can be violated at 
the local scale. General conditions are derived when the 
efficiency is higher than the noninteracting model. We 
also present an upper bound for efficiency which is lower 
than the Carnot bound. The proof is sketched in the Ap- 
pendix. In IIIC, we interpret some nontrivial features of 
the engine operation in terms of local spin temperatures. 
The final section IV summarises our findings. 



II. THE COUPLED QHE 

The working medium for our QHE consists of two spin- 
1/2 particles within the ID isotropic Heisenberg model 



2 



O 7 I O D 

Z.I -\- ZhS 


1 1 1 \ 






2J 


\1p+) = 

: |tt/ 




- 101)) 


2J-2B 


|00) 






-6J 






-|01)) 



FIG. 1. Energy eigenvalues and eigenstates of two-spins en- 
tangled model system 

[H, HH ■ The Hamiltonian is given by 

H = J(<7 1 .ct 2 + <j 2 .<j 1 ) + B(<t 1 z +ct 2 z ), (1) 

where er^ 2 ' = (<Ts , <7y , <ri^) are the Pauli matri- 
ces, J = J x — J y = J z is the exchange constant and 
B is the magnetic field along z-axis. Cases J > and 
J < correspond to antiferromagnetic and ferromagnetic 
interactions, respectively. In this paper, we consider an- 
tiferromagnetic case only. The energy eigenvalues of H 
are -6J, (2 J - 25), 2J and (2J + 25). If |0) and |1) 
represent the state of the spin along and opposite to the 
direction of the magnetic field respectively, then in the 
natural basis {|H), |10), |01), |00)}, we can write the 
density matrix as 

/ 7=P 1 |V^)(^-|+ J P2|00)(00|+P 3 |^ + )(^ + |+P4|ll)(H|, 

(2) 

where \ip±) = (| 10) ± |01))/V2 are the maximally en- 
tangled Bell states. The occupation probabilities of the 
system in the thermal state at temperature T are given 
by 
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where, Z = (1 + e 8J / T + e 2B / T + er 2B l T ) is the normal- 
isation constant. 



A. The heat cycle 

The four stages invloved in our quantum Otto cycle 
are described below: 

Stage 1 : the system with the external magnetic field at 



B\ attains thermal equilibrium with a bath of tempera- 
ture T\. Let occupation probabilities be pi, p%, P3, and 
Pi as tabulated above with T = T% and B = B\. Stage 2: 
the system is isolated from the hot bath and the magnetic 
field is changed from B\ to B 2 by an adiabatic process. 
According to quantum adiabatic theorem, the process 
should be slow enough to maintain the individual occu- 
pation probability of each energy level. Stage 3: the 
system is brought in thermal contact with a cold bath 
at temperature Ti. Upon attaining equilibrium with the 
bath, the occupation probabilities become p[, p' 2 , p' 3 , and 
P4 corresponding to the thermal state with T = T% and 
B = B 2 . On the average, the system gives off heat to the 
bath. Stage 4 ■ the system is removed from the cold bath 
and undergoes another quantum adiabatic process which 
changes the magnetic field from B 2 to B\ but keeps the 
probabilities p[, p' 2l p' 3 , and p' 4 unaffected. Finally, the 
system is brought back to touch the hot bath. On the 
average, heat is absorbed from the bath and the system 
returns to its initial state. 

The heat transfered in Stage 1 and in Stage 3 of the 
cycle respectively is 

Qi = £>iOi-p0 (7) 

i 

= 8 J(pi - Pi) + 2B 1 {p' 2 -p 2 + P4 ~ ti), (8) 

and 

Q 2 = J2E'M~Pi) (9) 

i 

= -8J(pJ -pi) - 2B 2 (p' 2 -p 2 + Pi -p , i ). (10) 

In the above, Ei and E[ (i — 1,2,3,4) are the energy 
eigenvalues of the system in Stage 1 and Stage 3 respec- 
tively. Qi > and Q 2 < corresponds to absorption 
of heat from hot bath and release of heat to cold bath 
respectively. Comparing the equations for heat trans- 
fer between the system and the reservoirs, Eqs. (J5J) and 
([TO]) , the quantity of heat 8J(p[ — pi) appears in both 
the equations. Obviously, this term is absent in the un- 
coupled case for which J — 0. As will be shown below, 
the sign (±) of this term determines whether the effi- 
ciency in the coupled case will be higher or lower than 
the uncoupled case. 

The work is done in Stage 2 and Stage 4 when the en- 
ergy levels are changed at fixed occupation probabilities. 
The net work done by the QHE is 

W = Q 1 + Q 2 = 2(5! - B 2 )(p' 2 -P2+P4- Pi). (11) 

Note that W > corresponds to work performed by the 
system. 

III. THE LOCAL DESCRIPTION 

In this section, we discuss how the individual spins in 
the system undergo the cycle. Again, let Q\ 2 and g' 12 rep- 
resent the thermal states in the natural basis when the 
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system is in equilibrium in Stage 1 and Stage 3 respec- 
tively. Explicitly, the density matrices are 
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(13) 



Let Qi and g 2 be the reduced density matrices in Stage 1 
for the first and the second spin, respectively. Then from 
the normalization constraints, Yijpi — = 1, we get 



Qi = Q2 



(P2-P4) 
2 







(P2-P4) 



(14) 



Similarly in Stage 3, the reduced density matrices for the 
first and second spin are 



Qi = Q 2 



2 









1 I (P2-P4) 

2 T 2 



(15) 



Since the applied magnetic field is the same for each spin, 
their local Hamiltonian is also same. Let Hi and H! be 
the local Hamiltonians for individual spins with eigen- 
values £?i) and (B 2 ,—B 2 ) in Stage 1 and Stage 3 
respectively. The heat transferred locally between one 
spin and a reservoir is given by 

<l\ = B l(p' 2 ~P2 +Pi -p'a), (16) 
q 2 = -B 2 (p 2 - P2+P4- P4), (17) 

for the hot and the cold reservoir, respectively. So we get 
the net work done by an individual spin as 



w = qi + q 2 = {Bi - B 2 ){p' 2 - 
From Eqs. (QJ§ and JTTJ 

W = 2w. 



P2+Pi-p' 4 )- (18) 



(19) 



Thus the total work performed is the sum of work ob- 
tained from the two spins locally. 

Further, the total heat absorbed by the system can be 
written as 

Qi = 8J(pi-pi)+2gi, (20) 

and similarly for the heat released to the cold bath is 

Q2 = -8J{p' 1 ~ Pl ) + 2q 2 . (21) 

Now it can be seen that because the work is done only due 
to change in local hamiltonians, so only the part of the 
heat which is absorbed locally by a spin can be converted 



into heat. The part 8J(p[ — p\) cannot potentially be 
converted into work due to the nature of the adiabtic 
process involved and is transferred directly between the 
reservoirs. But it may not be transfered only from the 
hot to the cold bath, in which case it may be regarded 
like a heat leakage term. In fact, the flow of this heat can 
be in the opposite direction which is directly related to 
the enhancement of efficiency due to coupling, as shown 
below. 

In the following, we consider two cases whereby mag- 
netic field may be decreased or alternately, increased in 
Stage 2. It will be seen that the second case is feasible 
only in the presence of interactions, J ^ 0. In the first 
case when J = 0, the above equations go back to those 
for Kieu's model with two uncoupled spins where an en- 
gine operation is obtained given T\ > T 2 and B\ > B 2 
with the additional condition B 2 /T 2 > B\JT\. 



A. The case Bi > B 2 

From Eq. (| 1 1 1) , the condition that the work performed 
be positive (W > 0) is given by 



(P2 - Pi) > (P2 - Pi)- 



(22) 



Secondly, for the heat to be absorbed from the hot bath 
(Qi > 0), from Eq. ((8]) we have one of the following two 
possibilities: (i) p[ > p\ or (ii) p[ < p±. Alongwith the 
possibility (ii), we must also have (p' 2 — p 2 + P4 — p' 4 ) > 
(4J/B 1 )(p 1 -pi). Now we rewrite Eq. flHJ as 



Qi = 8J(p[-pi) 



WBi 



{B x -B 2 ) 



(23) 



or 8J(p' 1 — pi) — Qi(l — 77/770)7 where 77 = W/Qi is the 
efficiency of the coupled engine and 770 = {B\ — B 2 )jB\ 
is the efficiency of the uncoupled i.e. J = case. Thus 
for J > 0, if p[ > pi, then 77 < 770, or the presence of 
coupling between the spins decreases the efficiency from 
its value 770. The global efficiency is equal to the local 
efficiency in two situations, when J = or p\ = p\. 

On the other hand, if p[ < p\. then it is possible that 
the efficiency of the coupled engine can be higher than 
the uncoupled case. Using the latter condition with Eq. 
|), we have 



(Pg -Pi) > (P2 ~ Pi) 



Pi 



Pi 



(24) 



From the explicit expressions for the probabilities, the 
above inequality can be simplified to give 



B2 Bi 
T 2 TV 



(25) 



Thus we see that the above condition which is necessary 
to extract work in the J = model is also the condition 
for the coupled case to obtain an efficiency higher than 
770. But additionally, for a set of given values of T\, T 2 , 
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FIG. 2. Efficiency versus the coupling constant J, for Bi > B2 
case, for values Si = 4, B 2 = 3, Ti = 1 and T 2 = 0.5. The 
uncoupled case corresponds to 770 = 1 — S2/B1 — 0.25. The 
upper curve denotes the bound for efficiency from Eq. (|26|) . 
Beyond the maximum shown value of J, the efficiency drops 
monotonically. The Carnot limit is 0.5. 



B\ and B 2 , there is a maximum value of J beyond which 
the efficiency drops below the 770 value. See Fig. [2J 

The reason for the lowering of efficiency when p\ < p\ , 
is that the term 8J(p' 1 — pt) is positive and it acts like 
heat leakage term which reduces the efficiency. On the 
other hand, when p\ > p[, this term is negative which 
means that although each spin locally absorbs heat equal 
to qi from the hot bath, due to interaction the effective 
total heat absorbed by the two-spin system is less than 
2gi, which raises the efficiency for a given quantity of 
the work performed. It is interesting to know how much 
maximum gain in efficiency is possible for a given set of 
parameters . We have proved an upper bound for the 
global efficiency, given by 



1 - Ba/Bi 
1 - U/Bi 



(26) 



where r\ c = I — T2/T1 is the Carnot bound. Also for 
77 > 770, we have the condition B\ > 4 J. This implies that 
the ordering of energy levels which gives an enhancement 
of efficiency (over the uncoupled model) is: 



(2 J - 2B t ) < — 6J < 2 J < (2 J + 2Bi) 



(27) 



and which after the first quantum adiabatic process, be- 
comes 

(2J-2B 2 ) < -6J<2J< (2J + 2B 2 ). (28) 
The proof of Eq. (f2"6"|) is given in the Appendix. 

B. THE CASE B 2 > Si 

In this case, during the first quantum adiabatic pro- 
cess, the magnetic field is increased from its value Bi to 
B 2 . If there is no interaction between the spins, the sys- 
tem cannot work as an engine in this case because the 



condition W > will not be satisfied [|| . The conditions 
Ti > T 2 and B 2 > Bi directly lead to 



Pi > Pi, 

P3 > P3- 



(29) 
(30) 



Further, the positive work condition implies (p 2 — p'4) < 
(p 2 — p^, which alongwith (|2"9")l gives 



P2 > P 2 - 



(31) 



The normalisation of probabilities and the above three 
conditions Eqs. (|29]). ([30]) and ([31]) together give 



Pi > pi- 



rn) 



These are the necessary conditions for the system to work 
as an engine given that Ti > T 2 and B 2 > B\ . According 
to Eq. (jTSJ) , the local work should be positive. This 
yields qi < and q 2 > 0. This means locally the heat is 
absorbed from the cold bath and given to the hot bath. 
Also the local efficiency is 



W 

<T2 



1 



Bi 



(33) 



Thus locally, the spins operate counter to the global tem- 
perature gradient present due to T\ > T 2 . But globally 
we do have Q\ > and Q 2 < 0. Thus the function of the 
two-spin engine is consistent with the second law of ther- 
modynamics, although locally we seem to have a viola- 
tion of the same. This apparent contradiction is resolved 
below using the concept of local effective temperatures. 



C. Local temperatures 

Now each spin in the 2-spin system can be assigned 
a local effective temperature, corresponding to its local 
thermal state or the reduced density matrix [24- 26] . This 
is true regardless of the state of the total system. Partic- 
ularly, in stages 1 and 3 of the cycle, from Eqs. (|14]|15[) 
alongwith local Hamiltonian, we get the local tempera- 
tures as 



T[ = 2B X lo 



T' 2 = 2B 2 [ log 



(l+j>4 -P2) 



- I 



- I 



(34) 



(35) 



The important fact is that in the presence of interactions, 
the local temperatures are different from the correspond- 
ing bath temperatures. Thus T{ ^ T Y and V 2 ^ T 2 if 
J 0. Further, since the work in our heat cycle is done 
only locally, the total work by the system can be regarded 
as equal to the work by two independent spins operating 
between their effective temperatures. 
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(i) Engine working in B\ > B 2 : the positive work con- 
dition for a single spin is given by 



B 2 B\ 

rpl ^ rpl ■ 

1 2 1 1 



Since Bi > B 2 , we get 



T[ > T' 2 . 



(36) 



(37) 



and T 2 



To and we have the result 



At J = 0, T{ = T\ 
of Kieu's model [4]. 

(ii) Engine working in B 2 > B\\ in this model, the 
positive work condition is satisfied only when 

£ > (38) 

Thus in this case T 2 > T[. Moreover, it can be shown 
from the definitions (|34)) and f|35[) that for both the cases, 
T{ > Ti and T-^ > T 2 . Finally, based on local temper- 
atures, the counter-intuitive mechanism which leads in 
case (ii) to q\ < and q 2 > can be justified as follows. 
For B 2 > B\, due to the first adiabatic process, the local 
temperature increases from T[ to T[(B 2 / B\). After con- 
tact with the cold bath, the local temperature becomes 
T 2 , which due to condition (f38)) is more than T[{B 2 j B\). 
Thus heat should flow from the cold bath to the spin or 
q 2 > 0. Similar considerations lead to rejection of heat 
by the spin at the hot bath or q\ < 0. 



IV. SUMMARY 



APPENDIX 
Upper bound for global efficiency 

We consider the case of the engine working in the range 
B\ > B 2 . The condition to get a higher efficiency as 
compared to uncoupled model is the case (ii) discussed 
in Section IIIA and is given by 



Pi > Pi- 

From the condition B 2 /T 2 > Bi/T x (Eq. ([25 

P3 > P3, 
Pi > Pi- 

Then normalisation of the probabilities gives 
P2 > Pi- 

From Eqs. and (g2]), we have 

P2 



4>- 

Pi Pi 



(39) 
we get 

(40) 
(41) 

(42) 
(43) 



which simplifies to 

e (B 2 -4J)/T 2 > e (B 1 -4.J)/T K 



(44) 



Fig. [3] shows three possible ways of arranging the energy 
levels (2 J — 2£?i) and — 6 J relative to the level (2 J — 
2B 2 ) resulting from the first quantum adiabatic process. 
Equivalently, Eq. (|44|) is of the form e x > e v , which may 
be satisfied in one of the following three ways: 



A model of coupled spins is used as working medium 
to realise a quantum Otto engine. The conditions for the 
efficiency to be higher than the non-interacting case are 
found. The antiferromagnetic interaction between the 
spins allows a fraction of the total heat to flow from cold 
bath to hot bath provided the total heat should flow in 
in a direction suggested by global temperature gradient. 
This mechanism increases the efficiency of the system 
compared to non-interacting spins case. A tighter upper 
bound for the efficiency is found which is lower than the 
Carnot value. The system can also work as a heat engine 
even if it undergoes an adiabatic compression (B 2 > B{) 
in the second stage of the cycle. Here we have observed 
an interesting mode of operation using the reduced den- 
sity matrix whereby each spin absorbs heat from the cold 
bath and rejects some heat to the hot bath while perform- 
ing a net work. This feature is also confirmed from the 
analysis of local effective temperatures of the spins. 
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6 J 2J~2B 2 2J-2B 2 

2J-2B 1 6 J 

2J-2B 2 

2J-2Bi 6 J 2J-2B 1 

(a) (b) (c) 

FIG. 3. Three possible configurations of energy levels with 
eigenvalues — 6J, (2J — 2B\) and the level (2J — 2B2) result- 
ing from the first quantum adiabatic process whereby B\ is 
changed to a lower value B^- Only case (a) is possible as 
discussed in the Appendix. 

Case (a) represents y > 0, x > and so x > y. This 
implies, Bx > 4J and B 2 > 4 J. 

Case (b) represents x < 0, y < and |a;| < \y\. This 
implies Bi < 4 J, B 2 < 4 J, but due to the fact T 2 /Ti < 1, 
we obtain B\ < B 2 which leads to a contradiction. 

Case (c) represents y < and x > 0. This possibility 
is also similarly ruled out. 
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So the only possibility is case (a) representing the fact 
that the energy levels (2J—2Bx) and (2J—2B 2 ) lie below 
the level — 6 J when the coupled engine gives a higher 
efficiency than the uncoupled case. 

When the inequality (|44l) holds, we can write 



B 2 - 4 J Bi - 4 J 



To 



Zi 



(45) 



Since £?i > 4 J , B 2 > 4 J and Ti > T 2 , we get 



Vo , _ 1 _ ^2 
1 - 4J/Si ?/c ~ Ti ' 



(46) 



where ?y = l — B^/Bx. Now the global efficiency defined 



as 77 = W/Qx, can be written as 



jj( Pl -p' 1 ) 



(47) 



From the inequalities between the probabilities (Eqs. 
(|39 ]) ,(|4T j) and (|42]>). it follows that (p x - pi) < (p 4 - 
+P2 — Pi)- Therefore, we finally obtain that when the 
efficiency is higher than the uncoupled case (or the lower 
bound is 770)1 then an upper bound for efficiency is given 
by 



77 < 



'/o 



1 - AJ/Bx 



< t]c 



(48) 



When J — 0, we have 77 = 770. A similar kind of proof 
can be constructed for the case B 2 > Bx- Interestingly, 
the same bound as Eq. (|4*5|) is obtained. 
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